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ABSTRACT 

Transport equations which describe the flow of holes and electrons in the 
heavily doped regions of a solar cell are presented in a form that is suitable 
for device modeling. Two experimentally determinable parameters, the effective 
bandgap shrinkage and the effective asymmetry factor are required to completely 
model the cell in these regions. Nevertheless, a knowledge of only the 
effective bandgap shrinkage is sufficient to model the terminal characteristics 
of the cell. The results of computer simulations of the effects of heavy 
doping are presented. The insensitivity of the terminal characteristics to 
the choice of effective asymmetry factor is shown along with the sensitivity 
of the electric field and quas i -el ectri c fields to this parameter. The 
dependence of the terminal characteristics on the effective bandgap shrinkage 
is also presented. 

INTRODUCTION 


The need for accurate modeling of the physical phenomenon which control 
the behavior of solar cells has long been recognized as being necessary to 
the development of high-efficiency solar cells. This paper is concerned with 
the accurate modeling of heavy doping effects and the proper technique for 
incorporating experimentally determined heavy doping parameters into a 
numerical solar cell model. 

Since both the emitter and the high-low junction regions of a high- 
efficiency solar cell are heavily doped, the performance of these regions is 
controlled by the various heavy doping effects which occur. Among the con- 
siderations for correct modeling of these regions are: 1) the regions are 

degenerate and, therefore, require the use of Fermi -D i rac statistics, 2) the 
dens i ty-of-states function is expected to have a change of shape due to band 
tailing or the formation of an impurity band, 3) the original energy states 
are expected to be shifted in energy due to free carrier- ion ized impurity 
interactions, and correlation effects. Recent papers have presented 
theoretical calculations of the magnitude of these various phenomenon, 
references 1 to 3. The problem, however, is that the present state-of-the 
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theory is not sufficiently developed. For example, the shape of the density- 
of-states functions is unknown, and without this information one is unable 
to compute the parameters which appear in the formulation of the transport 
equations in the presence of heavy doping effects. 

NOTATION 


6 

AE, 


D ,D 
P n 

i—i * 

J ,J 
P n 


••to 

n ie 

P, n 

q 

T 

V 


p >p 

p ’ n 


n c 

v V ,v c 


Equi 1 ibri 


effective asymmetry factor 
bandgap 

bandgap shrinkage 
effective bandgap shrinkage 
hole (electron) diffusion coefficient 
hole (electron) current density 
Boltzmann's constant 

intrinsic carrier concentration in a semiconductor with an 
unperturbed band structure 

effective intrinsic carrier concentration 

hole (electron) carrier concentration 

magnitude of the electronic charge 

absolute temperature (Kelvin) 

electrostatic potential 

hole (electron) mobility 

electron affinity 

The gradient operator with q (n c ) held constant 
urn values of the parameters listed above are denoted by a superscript 
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TRANSPORT EQUATIONS 


If the existence of a mobility edge separating localized from delocalized 
states is postulated, transport equations for materials with a position- 
dependent band structure can be written. These equations in the form given 
by Marshak and van Vliet (refs. 4, 5) are: 
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are the dens i ty-of-states effects. 

Without a knowledge of the shape of the dens i ty-of-state function and 
the location of the band edges these equations cannot be utilized for device 
modeling. Others have formulated transport equations under more restrictive 
assumptions, such as the rigid band approximation (ref. 6) or under conditions 
for which Boltzmann statistics apply, (ref. 7). While these formulations 
allow one to perform numerical calculations, they are restricted to cases 
which are known to be violated in the emitter and the back surface field 
region of high efficiency solar cells. 
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A recent reformulation of the transport equation by Lundstrom, Schwartz 
and Gray (refs. 8, 9) provides a means of avoiding these complications while 
retaining a precise model of the effects of heavy doping on the transport 
equations, and consequently upon the behavior of a solar cell (refs. 8, 9). 
They have shown that equations 1 and 2 can be rewritten as 


- kT„ p V p , (7) 


k \,V < 8 > 

where 
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is the effective bandgap shrinkage and 

Ay + 0 
A n 


is the effective asymmetry factor. 


( 10 ) 


The parameters, 9 and 0 , which contain the information on the position 
dependence of the dens i ty-of-ltates and the influence of Fermi-Dirac statistics 
are defined by, 


and 



kT 







( 11 ) 


( 12 ) 


The significance of recasting- the Marshak and van VI iet equations in 
the form shown by equations (7) and {5) is that the parameters which appear in 
equations (7) and (8) , namely Aq and y, are experimentally determinable. These 
experimentally determinable parameters contain the information on the band 
structure and the influence of Fermi-Dirac statistics necessary to model 
heavily doped semiconductor devices. In addition, the simple, Bol tzmann-1 ike 
form of these equations facilitates their use in device modeling. 
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The effective gap shrinkage, Aq, is the parameter inferred from typical, 
electrical measurements of heavy doping effects (refs, 10 to 12). Although 
this parameter is frequently referred to as the bandgap shrinkage, equation 9 
shows that Aq also accounts for the modified dens i ty-of-states and the 
influence of Fermi -Dirac statistics. The effective gap shrinkage can be 
related to the effective intrinsic carrier concentration as (refs 8, 9). 


2 2 AG/kT 

n. = n. e 
le io 


(13) 


where n. is the intrinsic carrier concentration for a material with 
unperturbed band structure. 


The parameter y, which we refer to as the effective asymmetry factor, 
accounts for the asymmetry in the shift of the band edges as well as the 
modified dens i ty-of-states and the influence of Fermi -Dirac statistics. Like 
Aq, y can be determined from electrical measurements, but unlike Aq, no 
measurements of y have been reported. 

Equations (7) and (8) can be used by the semiconductor device modeler to 
account for the complicated effects associated with heavy impurity doping. 
Unfortunately, information on Aq is quite limited, and no information on y 
exists. Fortunately, the terminal characteristics of the device are 
insensitive to y when the heavily doped regions are in low injection, quasi- 
neutral, and the dopants fully ionized, (refs. 8, 9). While it is true that 
the terminal characteristics of the device are unaffected by the choice of y, 
it is not true for some of the internal parameters of the cell. For example, 
the built-in potential, electric field and quasi-electric fields are strongly 
affected by y. 

The remainder of this paper is devoted to showing the results of computer 
simulations of solar cells to illustrate the use of the heavy doping model 
described above. The sensitivity of solar cell performance to the two heavy 
doping parameters, y and Aq is considered. 


SENSITIVITY OF SOLAR CELL PERFORMANCE TO y 

A set of computer simulations was performed in order to test the 
sensitivity of the device model to the choice of the effective asymmetry 
factor. x For these simulations, a p + nn + solar cell, described by the parameters 
listed in Table I , was chosen. The doping profile of the p-n junction in 


The numerical model used for these simulations is described in reference 9. 
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this device is shown in figure 1. Complete ionization of the impurities was 
assumed in all simulations. Two different choices for y were considered. In 
Case I, we set y = 0.5, and in Case II we put all of the effective band edge 
shift into the majoirty carrier band (i.e. y = 0 for p-type, y = 1 for n-type) . 
Comparison of the computed, one sun IV curves, shown in Table II, shows that 
the terminal characteristics of the cell are relatively unaffected by the 
choice for y. Similar agreement was observed when the spectral responses for 
the two cases were compared.. 

The comparisons of the computed I -V curves for the two cases confirms 
the prediction (refs. 9, 10) that the terminal characteristics of typical 
heavily doped semiconductor devices are not sensitive to the choice for y. 
Figure 2, however, which compares the open-circuit electric fields within 
the p-n junction region for the two cases, shows that the electric field is 
quite sensitive to the choice of y. The quasi-electric fields for holes and 
electrons are shown in figure 3. The two choices of y are seen to result in 
very different quasi -el ectric fields, nevertheless, the effective fields which 
act on the carriers are not affected by y. The effective fields are the sums 
of the electric and quasi -electri c fields for each carrier. Figures 2 and 3 
show that the effective fields for holes and electrons are nearly identical 
for the two cases. 

The open-circuit carrier concentrations within the p-n junction region 
are plotted in figure 4. As expected, (refs. 8, 9) the carrier concentrations 
are not sensitive to the choice of y, Since both the effective fields and 
carrier concentrations are independent of y, the current densities are also 
expected to be unaffected by y. Figure 5 shows that the expected result is 
observed. 


SENSITIVITY OF SOLAR CELL PERFORMANCE TO A Q 

The solar cell simulations showed that the terminal IV characteristics 
of typical silicon solar cells are, as expected (refs. 9, 10) not sensitive 
to the choice of y. The performance , of silicon solar cells, however, is quite 
sensitive to the effective bandgap shrinkage, Aq. In order to study the 
effect of Aq, we compared the simulated performance of the p + nn + solar cell 
with the simulated performance of an identical cell in which bandgap 
narrowing was suppressed (i.e., Aq = 0) . For these comparisons we chose 
y = 0.5 (all the results to be discussed, however, are independent of the 
choice of y) . 

The computed current versus voltage curves for the p + nn + cell with and 
without bandgap narrowing are compared in figure 6 which shows the wel 1 -known 
(refs. 13, 14) reduction in open-circuit voltage, Vqq, caused by bandgap 
narrowing. The cause for the reduction in open-circuit voltage can be as- 
certained by examining figure 7, a plot of the open-circuit energy band 
diagram for the p-n junction region of the cell. This energy band diagram was 
computed by employing the rigid band approximation and is for illustrative 
purposes only. A reduction of the potential barrier for electrons injected 
into the heavily doped p-type region is shown by figure 7. This reduction of 
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the potential barrier is caused by Aq and it is not affected by y. The 
reduced potential barrier results in increased injection of minority carriers 
into the heavily doped regions which increases the recombination currents 
associated with these portions of the device and, therefore, reduces the 
open-circuit voltage. 

The computed characteristics of the p + nn + cell, with and without bandgap 
narrowing, are summarized in Table III. Because this cell was designed 
(ref. 15 ) to minimize the deliterious effects associated with heavy impurity 
doping, the open-circuit voltage of the cell is not severely degraded by 
bandgap narrowing. In addition, since the short-circuit current of the cell 
is not significantly affected by bandgap narrowing, the fill factor of the 
cell is actually improved by bandgap narrowing. This effect, clearly shown 
by figure 6, minimizes the degradation of cell efficiency caused by bandgap 
narrowing. 

SUMMARY AND CONCLUSIONS 


A set of transport equations for use in analyzing heavily doped semi- 
conductor devices has been considered in this paper. The equations were first 
presented, and the interpretation of the two parameters, Aq and y, used to 
describe the effects associated with heavy impurity doping was discussed. The 
influence that these two parameters have on solar cell performance was then 
illustrated by the results of computer simulation. 

The most important conclusion reached in this paper is that heavily 
doped semiconductor devices can be modeled in terms of two experimentally 
determinable parameters. The effects associated with heavy impurity doping 
(i.e., band edge shifts, changes in the dens i ty-of-states and the influence 
of degenerate statistics) can be modeled accurately by using the two parameters 
in simple, Bol tzmann-1 ike transport equations. To model heavily doped semi- 
conductor devices, there is no need to introduce additional, unproven 
assumptions as is done in the frequently used rigid band approximation. In 
addition, we have demonstrated that the terminal I -V characteristics of 
typical, high-efficiency, single-crystal, silicon solar cells are not 
sensitive to the effective asymmetry factor, y. 
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TABLE I PARAMETERS USED TO MODEL THE p + nn + CELL. 


Temperature 
Device thickness 
Base dop \ n g den s i ty 
PN junction depth 
High-Low junction depth 
P + surface concentration 
N + surface concentration 
Doping prof i les 
SRH lifetime parameter, t 
SRH lifetime parameter, x nQ 
Auger coefficient, A 
Auger coefficient, A^ 

Front surface recombination velocity 

Back contact 

Solar spectrum 

Optical reflection losses 

Bandgap narrowing model ( Aq ) 

Effective asymmetry (y) : 

CASE I 
CASE l\ 


27°C 
300 um 

5 x lo’^cm" 3 
0 . 35um 
1 .Oym 

20 -T 
1 .5 x 10 cm 5 

3.0 x 10 2 °cm' 3 

complementary error function 

400ysec 

400ysec 

o « m“32 6 -1 

9.9x10 cm sec 

2.8 x 10 ^cm^sec ^ 

5000 cm/sec 

ohmi c 

AMI 

0 

Slotboom and DeGraaf 

Y = 0.5 

y = 0 for p-type 
y = 1 for n-type 
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TABLE II. COMPARISON OF 1 SUN l-V CURVES FOR TWO EFFECTIVE 
ASYMMETRY FACTOR MODELS. 


Vol tage 

Current Density (MA/Cm 2 ) 

Deviation 


Case 1 

Case II 


0 

37-4273 

37.4245 

7.5 X lo" 5 

0.1 

37.4259 

37.4232 

ro 

X 

0 

1 

0.2 

37.4243 

37.4216 

7.2 X lo' 5 

0.3 

37.4219 

37.4192 

7.2 x 10 -5 

0.4 

37.3932 

37.919 

3.5 x 10" 5 

0.5 

36.3423 

36.3910 

-1 .3 x 10‘ 3 

0.6 

14.0828 

14.9401 

-6.1 x TO -2 

0.62 

6.1596 

4.5492 

2.6 x lo' 1 

^OC Case I 

= 0.61512 volts 



V 0C Case 1 

( = 0 .61632 vol ts 




TABLE III. COMPARISON OF p + nn + CELL PERFORMANCE WITH AND WITHOUT 
BANDGAP NARROWING. 


Short-Circuit Current 
Collection Efficiency 
Open-Circuit Voltage 
Fill Factor 
Eff i ciency 


With BGN 

37*43 mA/cm 2 

0.98 

0 .615 Volts 
0.80 
0.192 


Without BGN 

37.59 mA/cm 2 
0.99 

0.676 Volts 
0.76 
0.201 
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Electric Field (V/cm) Doping Density (per cm 



Position (urn) 

Figure 1. Doping profile for the p-n junction 
region of the p+nn+ solar cell. 



Position (wm) 

Figure 2. Open-circuit electric field within 
the p-n junction region of the 
p + nn + cel 1 . 
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Quasi-Electric Field (V/cm) 



Position (»m) 

Figure 3a. Quasi-electric fields in the p-n 
junction region of the p + nn + cell 
(Case 1 ) . 



Position (nm) 

Figure 3b. Quas i -el ect r i c field within the 
p-n junction region of the p + nn + 
cell (Case II). The electron 
quas i -el ect r i c field is zero 
for this case. 


48 




Current Density (A/cm x 10 ) Carrier Concentration (per cm 5 ) 


holes Case I 

• • • * holes Case II 

------- electrons Case I 

o o o o electrons Case II 


0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0 

Position (pm) 

Figure k. Open-circuit carrier concentrations 
within the p-n junction region of 
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- 0 . 20-1 . . . . . . . . ■ 1 - 

0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0 

Position (pm) 

Figure 5. Open-circuit current densities 

in the p-n junction region of the 
p + nn + cell. 
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Figure 6. Current-voltage characteristics 
for a p + nn + solar cell with and 
without bandgap narrowing. 
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Figure 7. Open-circuit energy band diagram 
for the p + nn + cell with bandgap 
narrowing. 




